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Abstract 

We study phase transitions in uniformly frustrated S'C/(A^)-symmetric (2 + e)- 
dimensional lattice models describing type-II superconductors near the upper 
critical magnetic field H C 2(T). The low-temperature renormalization-group 
approach is employed for calculating the beta-function j3(T, f) with / an 
arbitrary rational magnetic frustration. The phase-boundary line H C 2(T) is 
the ultraviolet-stable fixed point found from the equation f3(T, f) = 0, the 
corresponding critical exponents being identical to those of the non-frustrated 
continuum system. The critical properties of the S[/(iV)-symmetric complex 
Ginzburg-Landau (GL) model are then examined in (4 + e) dimensions. The 
possibility of a continuous phase transition into the mixed state in such a 
model is suggested. 

PACS numbers: 74.60Ge, 74.60Ec 



1 



Typeset using REVTeX 



I. INTRODUCTION 



The challenging problem of the critical behavior of a type-II superconductor near the 
upper critical magnetic field H c2 (T) has a rich and long history going back to the seminal 
paper by Eilenberger [|l|] (see also @). It has been recognized for some time that an external 
magnetic field drastically changes the critical properties of superconductors. The magnetic 
field hinders the growth of the thermal fluctuations in the plane perpendicular to H, since 



the growth of their correlation length is restricted by the magnetic length scale I = ^Jhc/eH 
which is much shorter than the coherence length £. This effect of dimensional reduction 
results in an enhancement of the longitudinal fluctuations leading, in particular, to the 
increase of the lower critical dimension from 2 to 4. 

If critical fluctuations are ignored, the uniform frustration (in the language of spin models 
of the vortex lattice) eventually leads to a continuous phase transition into the Abrikosov flux 
lattice state. On the one hand, in contrast to mean-field theory, the standard renormalization 
group (RG) approach in 6 — e dimensions, in fact, failed so far to yield insight on the nature 
of the phase transition due to the appearence of an infinite number of invariant charges 
(relevant scaling variables) inherent to the non-renormalizable scalar 4 field theory in a 
field [||]. Physical arguments (presented e.g. in 0) support the existence of a first-order 
melting transition for the flux lattice. On the other hand, however, the conventional l/N- 
expansion, when applied to the S?7(iV)-symmetric Ginzburg-Landau (GL) model with an 
iV-component order parameter, gives a second-order phase transition above four dimensions 
1^-0. The question of the nature of the phase transition from the normal into the mixed 
state of a type-II superconductor remains therefore an open problem, even before the effects 
of impurities and the topical question of the vortex glass are to be considered. In this paper 
we wish to investigate what is to some extent a new question, namely the nature of this 
phase transition when a periodic lattice potential is coupled to the superconducting order 
parameter. 

It was shown some time ago j7|J§ that the interaction between thermal fluctuations and 
the underlying crystal lattice can restore the phase transition into the mixed state. It is well 
known that the vortex lattice state has two broken symmetries: i) phase coherence, that 
is off-diagonal long-range order (ODLRO), and ii) translational symmetry. Including the 
interaction with an underlying lattice means that the translational symmetry is explicitly 
broken and lattice effects (for instance Harper's broadening and splitting of the Landau 
levels suppressing infinite degeneracy of the energy spectrum) become particular acute. 

Lattice models of superconductors were introduced long ago and are of great interest 
by their own accord, being employed to study phase transitions in uniformly frustrated 
XY models f9|, or in some artificial condensed matter structures like two-dimensional (2D) 
Josephson-j unction arrays subject to a magnetic field [10] or in opals filled with a semicon- 



ductor or a metal 11 



Naturally, according to a naive point of view, near criticality the correlation length £ 
diverges and the system should "forget" about the discretness of the underlying lattice; 
however, this turns out not to be the case. In fact, in the continuum limit some specific 
uniformly- frustrated lattice models have been shown not to be equivalent to the conven- 
tional GL model H. If lattice effects are taken into account, in the long- wavelength limit 
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one obtains an infinite set of effective (renormalizable) GL Hamiltonians corresponding to 
different rational values of the frustration / = $/$o> where $ = Ba 2 is the magnetic flux 
per plaquette and $ — hc/2e the elementary flux quantum ||[J. 

The main goal of the present paper is to consider the critical behavior of lattice models of 
type-II superconductors in the normal phase above H C 2(T). We shall apply the standard low- 
temperature RG approach that was first introduced and developed in fl2]]l3]| . The method 
was first exploited for studying critical phenomena in superconductors (with B — 0) in 



14| . We will show how the presence of a lattice induces a second-order phase transition in 
dimensions d = 2 + e, by giving explicitly the phase boundary in the low-temperature limit. 

The remainder of this paper is organized as follows. In Section II the critical behavior of 
the (2 + e)-dimensional SU (iV)-symmetric lattice model is considered by means of the low- 
temperature RG approach. Section III deals with the standard GL model in 4+e dimensions, 
therefore in the absence of the underlying lattice, and we show that the phase transition 
appears to be again of second order, although it remains impossible to reach a conclusion 
about the nature of the phase transition in d = 3. Section IV contains a discussion of our 
results and of related issues in which lattice effects may play an important role and some 
concluding remarks. 



II. VORTEX-LATTICE PHASE TRANSITION IN 2 + e DIMENSIONS 

A. Effective action and low-temperature expansion 

We begin by considering the classical Hamiltonian of the S"?7(V)-symmetric nonlinear 
sigma model, defined on a square lattice with periodic boundary conditions: 

H = -J \S?e W {iAij) ~ S;\ 2 (2.1) 

<i,j> 

This is indeed the lattice version of the GL Hamiltonian, Sf being a iV-component complex 
unit vector with the constraint S^S" 1 * — 1, a — 1, 2, N and J is a coupling constant. Here 
< ■ • • > indicates that the summation is over all nearest-neighboring sites, as usual. Aij is 
a bond angle such that the sum around a plaquette is given by 

E Aj=2nf (2.2) 

plaquette 

with f = p/q the so-called frustration; p and q are here mutually prime integers and is 
defined by 

2tt p 

= — — / dx^A^ (2-3) 

SPq Ji 

where A^ = (—By, 0, 0) is the vector potential of the uniform magnetic field B along the z 
axis. From eq.( |2.2| ) it follows that the magnetic flux through a plaquette is assumed to be 
a rational fraction of the magnetic flux quantum $ - 

Consider the partition function associated with eq. Q2.lD , namely 
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J RdPSidTVexpi-jrWStS? - 1) (2.4) 



Notice that the Hamiltonian in eq. (|2.1| ) has both a global nonabelian SU(N) symmetry and 
a local i7(l) one whilst the fixed-length constraint imposed on the local spins is 0(2N)- 
symmetric. 

The main steps in our calculations are as follows. To carry out the weak-coupling ex- 
pansion for the theory given by eq.( |2.1| ) a formal procedure based on integrating out the 
high-frequency components of the local spins Sf will be applied. We shall make use of the 



following parametrization of spin variables [15] 



rrfc 


= A 




h „b* 


4 


= ^k^k 


s? 





1,2 N-l 



ix\ exp(i(j) k ) (2.5) 

where ir h k are small and slowly- varying fluctuations about the iVth-component's direction. 
Substituting the representation eq. Q2.5Q into eq. (|2.4|) and integrating out the modulus \S k \, 
one arrives at the Lagrangian 

Z = j l[d N - 1 TT l d N - 1 TT*ex P {-H/T} 

i 

H = H + H int 

H = -J E kfexp^O-^r 

<i,j> 



H in t = -J E I V 1 " *i eX P(^ - & + A v) ~ V 1 " ^ ( 2 - 6 ) 

<i,j> 

Here the Jacobian factor equals unity. The advantages of this parametrization are quite 
evident because if T — > the spin-wave fields v:\ can be treated as free lattice fields. As 
these fields fluctuate, a contribution to the effective action for the phases 0j will arise. At 
the same time, the field <fi(i) does not react back onto 7if in the renormalization procedure. 

The calculations in the one-loop approximation may be readily carried out by means 
of the standard momentum-shell recursion relation technique developed in [16|] (see, also 



15,17]). To produce a systematic low-temperature perturbation theory one has to expand 



nonlinearities such as y 1 — nf in Hi nt in powers of 11? and to integrate out short- wavelength 
degrees of freedom. Let us decompose the Fourier-transformed spin field according to 

n b (q) = n b < (q) < q < A' 

vr 6 (g) = 7rt(g) A' < q < A (2.7) 

with the purpose of integrating out the short-wavelength fields vr^(g). Here A and A' are 
momentum cutoffs and q stands here for the reciprocal lattice coordinates. Calculations in 
the lowest order of perturbation theory yield 
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Hint _ 2 \ „„„f ± j. i A \fi 1 ^ „2 ^ 1 .- _2 



- £ 008(0, - & + ^)(1 - - < TT? > - - < 7T| >) (2.8) 

<ij> 



where the angular brackets < • • • > stand for the straight Gaussian integration over the 
modes ir>(q), and J is absorbed in T. From eq. ([2.8|) the temperature renormalization is 
easily seen to be 

Y, = ^(!- < < >) ( 2 - 9 ) 

Setting the magnetic field (or, equivalently, the frustration) to zero one may easily take the 
long-distance limit of the model and find the bare propagator for the 7r> (q) fields, namely 

G ab (q) =< 7r£(g)7r£(-g) >= ^5 ab (2.10) 

q A 

The most divergent part of the Goldstone propagator behaves like 1/q 2 for small q, giving 
rise to the dominant logarithmic term in < 7if >= (N — 1) hip-. Carrying out the trivial 
calculations as described in detail in [15|, one is lead to the familiar expression for the 
one-loop beta function 

/3(T) = (d-2)T- ( iV ~ 1 ) T2 (2.H) 

7T 



The problem now is to extend this approach to the uniformly frustrated model of eq. (|2.lQ 
that is when the magnetic field is switched back on. 



B. Lattice frustrated models 

Since we are, in this paper, interested in the critical behavior of the lattice frustrated 
model, we shall first of all note some properties of the Azbel-Harper-Hofstadter operator L 
|T8VpT] describing a 2D Bloch charged particle subject to a uniform magnetic field. Here L 
is the inverse of the quadratic part of the action Ho which one may formaly regard as the 
electron hopping Hamiltonian. Using a compact notation, the Hamiltonian Hq in eq.( |2.6|) 
can be rewritten as 

ffo = £**^» (2.12) 

n 

where by definition (due to the choice of the most conventional Landau gauge) L acts on 
the Bloch wave function in the following way 

Ln b (n) = exp(ik x )7rl(n — 1) + exp(— ik x )T[ b k {n + 1) 
+ 2cos(k y + 2n—)n b k (n) 

n b (n) = exp(ikn)n b k (n) (2.13) 

where n = 1, q is a coordinate in the magnetic cell and the wave vector k ranges over the 
reduced Brillouin zone: — n/a < k x < it /a, —n/qa < k y < n/qa. It is worth noting that if n 
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labels the components of a g-component "vector" 7r|(n), then the operator L acts as a g x g 
Hermitean matrix. In the case of g = 4 one has, for example 



( 2 cos(k x ) 
exp(ik y ) 


V exp(— ik y 



exp(-iky) 

2 cos(k x + 27T/") 

exp(ik y ) 









2008(^ + 4^/) 
exp(ifc 3/ ) 



exp(ifc ?/ ) ^ 


exp(— iky) 
2 cos{k x + 6nf) j 



(2.14) 



Therefore, this Hamiltonian Ho can be thought of as that corresponding to a particle hopping 
along q sites around a ring P2,25|. The spectrum of Hq is known to posess a finite set of 
q magnetic sub-bands, each state in these sub-bands being g-fold degenerate instead of the 



infinite degenerancy inherent to the continuous problem [23 



The magnetic translational symmetry properties result from the local gauge invariance of 
the lattice theory, eq. (|2.1| ). This implies that under the translation r — > r + a the propagator 
in eq.( p,10| ) transforms as follows 



Z7T 



G ab {r + a, r' + a) = exp{ — [B x a](r - r')}G ab {r, r') 

<P 



(2.15) 



From eq. (|2.15|) it follows that G a b(r, r') at coinciding points does not depend on r. On the 
other hand, for the propagator from eq.( [2.15|) one has the resolvent spectral decomposition 



G ab (v, r') = cU E E E ^nk*nka(r)*; ka (r') 

n=l a=l k 



(2.16) 



where n = l,.--,q labels the magnetic band number, and \I/nka(r) are eigenf unctions of 
L defined in the magnetic Brillouin zone, — it /a < k x < it /a and —irjqa < k y < n/qa, 
forming the basis of the g-dimensional irreducible projective representation of the magnetic 
translation group. The exact energy spectrum E nk does not depend on the quantum number 



a: 



q. 



There exists a remarkable property of the energy spectrum resulting from the local gauge 
invariance of the lattice theory, eq.( [2.1| ). From this local gauge invariance it follows that 
Enk depends on k only through the parameter A = 2 — cos(qk x ) — cos(qk y ) and of course on 

E n]l = E n {A,q) (2.17) 

From the point of view of the critical properties of the nonlinear cx-model under consideration, 
one can see that only states lying near the bottom of the lowest magnetic sub-band (n = 1) 
are relevant. Here it works in a similar way, as in the well-known lowest Landau level (LLL) 
projection approximation for the continuous theory. 

Within the effective mass approximation the energy spectrum as a function of a quasi- 
momentum near the band's bottom (A = 0) reads 



E 



ik 



eo(?) 



{{qk x f + {qkyf} 



(2.18) 
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with 60,2(9) being coefficients of the Taylor expansion of Ei^ in powers of k. Here e 2 {q)~ 1 is 
proportional to an effective mass in the (x, y)-plane, eo(q) being proportional to a chemical 
potential. It is here rather essential that the effective mass and chemical potential are some 
functions of q. 

In order to find j3(T) we have to calculate < 7if >, which is given by 



a=l 

= ? EE 3* ( 2 - 19 ) 



n=l k 



The factor q appears as a result of the degeneracy of the spectrum and summation over 
a. Making use the effective mass approximation for the lowest magnetic sub-band, and 



evaluating the integral over k in eq.( p.l9| ), we arrive at the following expression 

262(g)- 1 



(qk x ) 2 + (qk y ) 2 + ml 



2(N-l)q A 

= 1 , , - In 77- 2.20 

Here, Wq = 2e (q)/e 2 (q) denotes some effective mass. With the help of eq.( |2.9|) and of 
eq.( |2.20D , one can readily derive the beta-function in the one- loop approximation 



m Q) = (d- 2)T - 2{N ] ] f 2 ( 2 - 21 ) 

ne 2 {q) 



This expression is the lattice version of eq. (|2.11 ). It leads to the nontrivial ultraviolet-stable 
fixed point located at 



We see that, in fact, eq.( 2.22j ) gives the phase transition line H c2 (T) of the lattice model, 



eq.( J2.1|j . The peculiarity of the beta-function, eq.( |2.21p , lies in its explicit dependence on 



the frustration q revealing a multifractal type structure. This follows from the properties of 
the spectrum of a single-electron in a magnetic field and a periodic potential (leading to the 
so-called "Hofstadter's butterfly" multifractal structure [|T9| , pCi|| ) . From the physical point of 
view this peculiarity looks quite natural since the phase transition temperature must depend 
on the applied magnetic field. After setting q — 1, or H — 0, e 2 {q) becomes e2(l) = 2 and 
the function (3{T) coincides with eq.( [2~9"|) . 

The critical exponent for the correlation length v follows immediately by differentiating 
eq. (|2.21|) at the fixed point. Keeping only the leading term in e one obtains 



" = ~wh = 7 (2 ' 23) 
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To determine Fisher's critical exponent 77 one may apply the conventional momentum-shell 



RG approach to first order in e as described in detail in [16| , [17| . After carrying out the 



standard RG computational procedure, we are led to the following expression 

n{ T) = 2- d+ V N - 1 f (2.24) 
7re 2 (g) 



Note that eq.( ggg|) and eq.( ggg) yield the usual one- loop result r/ — 2 (n-i) ■ 

We have seen above that the (2 + e)-dimensional frustrated lattice model undergoes a 
continuous phase transition on the coexistance curve H c2 (T) and is characterized by universal 
critical exponents which we have calculated. A physical interpretation of this conclusion is 
quite simple: the infinite degenerancy inherent to a charged particle moving in a uniform 
magnetic field is lifted by a commensurate potential suppressing the dimensional reduction 
effect and favouring the flux lattice state ||. At least at the lowest order in e the frustration q 
drops out of the critical exponents, these being identical to those for the conventional 0(2A)- 
symmetric Heisenberg ferromagnet. For N = 1 the expressions obtained above obviously 
show a singularity in the factor j^zf, which reflects the special properties of the neutral 2D 
XY model. 



III. WEAK-COUPLING EXPANSION FOR THE GINZBURG-LANDAU MODEL 

The approach which has been developed until now may be extended to the continuous 
GL model subject to a uniform magnetic field B in d = 4 + e dimensions, in the spirit of the 



work of Lawrie and Athorne JT4],f^J • This model is described by the Hamiltonian 

/2tt 
d d x\(d, + i—A^ a \ 2 (3.1) 
<p 

where the summations over fi = l,...,d and a = 1,..,N are understood in eq. (|3.1| ). Here 
\I/ = \I/jv] is an A-component complex order parameter, the fixed-length constraint 

|\I/| 2 = 1 being imposed on the local fields. The vector potential is now taken within the 
symmetric gauge 

A = x r (3.2) 
where B is taken along the z axis. The partition function associated with eq.( |3.1| ) reads 

N 

D^ a D^* a exp{- 

0=1 

We shall make use of the same parametrization of \l/ a exploited in Section II 



r N TT 

Z = JU D^ a D^ a ex V (- Y )5m 2 - 1) (3.3) 



# 6 (r)=7T 6 (r) 6= 1,2,..., A- 1 
7r 2 (r) = 7T 6 (r)<(r) 



^(r) = Jl- vr 2 (r) exp(z0(r)) (3.4) 
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Inserting eq.( |3.4j ) into eq. ( |3.3|) and integrating out the \^n\, we obtain the effective Hamil- 
tonian convenient for generating the low-temperature expansion 



r 27T ?7T , 

H= d d x[\{d, + i—A,)Tt b \ 2 + {d^+—A,Y + {d,VT^ 1 ) 2 ] (3.5) 

J <P <&0 



Observe that no Jacobian arises from the elimination of |^jv|. To determine the temperature 
renormalization in the lowest order of perturbation theory, we have to take into account only 
the first two relevant terms in eq. ( |3.4j ) and disregard the third one describing the interaction 
between the transverse degrees of freedom 7T&. The RG equation obtained in this way is 

^ = ^(1- < ^(r) >.) (3.6) 

where < ■ • • > s stands for integrating over the short-wavelength transverse fields 7Tb. 
Notice that the non gauge-invariant correlation function 

G a6 (r,r / )=<7r a (r)7r,(r / )> (3.7) 
is the Green function of the d- dimensional Schrodinger operator 

(-^ - :?r AJ'GWr, r') = 5 ab 5(r - r') (3.8) 

The exact solution of eq. fl3.8|) in an arbitrary gauge reads |26" | 



2-d 



G ab (r,r') = 5 ab exp I -i— J dx^A^ J (Atx)~ du 



! : - z') 2 , I 



2 sinh(wa;/2) 



x exp{- v ~ 4 J 7 - | coth[^((a; - x') 2 + (y - y 1 ) 2 )]} (3.9) 

where to = is the cyclotron frequency (here we have set % = 1 and 2m = 1) and z and z' 
are (d — 2) -dimensional longitudinal coordinates. The integral in eq. (|3.9|) is taken over the 
straight line connecting the points r and r'. 

Near the phase boundary H c2 (T), only the lowest Landau level gives a dominant con- 
tribution to the gauge- invariant quantity G(r, r) entering eq. Q3.6p . The LLL approximation 
leads to a nice simplification of eq. (|3.9|) 



G ab {r, r') = 5 ab exp{-i^- f dx^ -%[{x- x') 2 + (y - y') 2 }} 

/. x d f°° , 2 ~ d r uuj iz — z') 2 , 

x 4vr "2 / duu— exp{ ^ ! -\ 3.10 

Jo 2 4m 

In going from eq.( [3.9|) to eq . ( p . 1 0|) we have made the following substitution in eq. (|3~^ ): 
coth(Mdj/2) — > 1 and sinh(-uc<j/2) — > |exp(wu;/2) (justified by the limit B — > oo of the LLL 
approximation) . 

To make further calculations more transparent, we will make use of the mixed coordinate- 
momentum representation for G ab (r,r'). After carrying out the Fourier transformation in 
the longitudinal variables z, z' eq. Q3.10| ) becomes 
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G ab (x,x',y,y',k) = 5 ab exp{-i^- dx^A^ - - x') 2 + (y - y') 2 ]} fc2 ^ ^ 2 (3.11) 



where k is a d — 2-dimensional vector and £ is a correlation length in the longitudinal 
directions. The important conclusion one may draw from eq. ( |3.11|) is that a dimensional 



reduction effect takes place in the physics of our model, eq. (|3.1| ). We have also set an 
irrelevant factor to unity. 



With the result of eq. (|3.11|) at hand, we can now readily evaluate < n (r) > s 



(N -l)u , A 

< 7r (r) > s = ^ln- (3.12) 

7T A 

Eq.( |3.6| ) and eq.( |3.12| ) yield the beta-function in the one- loop approximation as well as the 
phase boundary, namely 

P(T, B) = (d- 4)T - (iV ~ lV T 2 

7T 

It is remarkable that the calculation of the critical exponents, like the one carried out 
in Section II, leads now to the same universal values as given by eq. ( p.22j ) and eq. ( [2.24j ), 
independently of B. 

We have seen that in contrast to the results obtained e.g. in || for the (6 — e) -dimensional 
case, a continuous phase transition occurs in the GL model in 4 + e dimensions. It is 
described by the critical exponents of the 0(2iV)-symmetric Heisenberg ferromagnet (in the 
one-loop approximation). From the field-theoretical point of view we are now dealing with 
the dimensional reduction effect within the Abelian Higgs model defined by an Euclidean 
Lagrangian, eq. (|3.1| ) (or, equivalently, in the A^-component scalar QED), in a large external 
magnetic field. The analogous effect was recently shown to occur also in the conventional 
4D spinor QED subject to an extremely high magnetic field like in a vicinity of a neutron 
star 



IV. CONCLUSIONS 

It has been shown that the (2 + e) -dimensional 5'L r (A r )-symmetric uniformly-frustrated 
lattice spin model undergoes a second-order phase transition, described by the universal 
critical exponents of the 0(2AQ-symmetric Heisenberg ferromagnet irrespectively of the value 
of the frustration. Our calculations provide no evidence of a first-order phase transition. This 
result was found to hold both for the (2 + e) -dimensional lattice model as well as for the 
(4 + e) -dimensional continuos GL model in an applied magnetic field. The occurence of the 
second-order phase transition is in this case in contrast to the RG calculations carried out 
in (6 — e)-dimensions @. Contrary to the conventional RG approach based on the 4 -theory 
||, the low-temperature RG approach employed in this paper works for an arbitrary rational 
value of the frustration / = $/$o- R allows us to compute the phase transition boundary 
line H c2 (T), which is a very difficult calculation problem within the standard 4 -theory. In 
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the case of the lattice model the coexistence curve exhibits a very complicated structure 
(in fact, reminiscent of the "devil's staircase" || inherited from the Hofstadter's butterfly's 



structure p9|j20[|) which can be determined only numerically for all integer values of q. 



Our treatment can be easily extended to the large- N limit. But the situation here is 
somewhat unclear since different groups of researchers have come to different conclusions 
concerning the nature of the transition. On the one hand, our results - being in good 
agreement with those obtained within the -^-expansion by Radzihovsky |3j and by the 
authors of - give evidence in favour of a continuous phase transition. In particular, the 
non-trivial fixed point responsible for the second-order phase transition was found. Moreover, 
from the equation for the correlation length £ it follows that there is a divergence of £ at 
crticality. One the other hand, the calculation of the effective potential within the ~- 
expansion as presented in [ffl , gives evidence for a first-order phase transition. The origin of 
this discrepancy in conclusions is rather subtle and probably lies in non-commutativity of 
the large- N limit and the thermodynamic limit (see, for a discussion, ref. @). 

The main problem still open for future research is how to extend the results obtained by 
means of the low-temperature RG approach to the N — 1 case. Thus, we have seen that 
the underlying lattice plays an essential important role near the upper critical magnetic 
field H C 2(T) restoring (like weak disorder) the continuos phase transition into the flux- 
lattice state. This is indeed the main result of our paper and we have commented in the 
Introduction on possible direct applications of the model we have studied. 
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